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Abstract 

We compare form factors in sine-Gordon theory, obtained via the bootstrap, to 
finite volume matrix elements computed using the truncated conformal space ap- 
QJ I proach. For breather form factors, this is essentially a straightforward application of 



(N 
> 



a previously developed formalism that describes the volume dependence of operator 
matrix elements up to corrections exponentially decaying with the volume. In the 
case of solitons, it is necessary to generalize the formalism to include effects of non- 
diagonal scattering. In some cases it is also necessary to take into account some of 
■ the exponential corrections (so-called ^-terms) to get agreement with the numerical 

data. For almost all matrix elements the comparison is a success, with the puzzling 
^ exception of some breather matrix elements that contain disconnected pieces. We 

' also give a short discussion of the implications of the observed behaviour of /_f-terms 

on the determination of operator matrix elements from finite volume data, as occurs 
e.g. in the context of lattice field theory. 

> . 

H 1 Introduction 

, 

The matrix elements of local operators (form factors) are central objects in quantum field 
theory. In two-dimensional integrable quantum field theory the S matrix can be obtained 
exactly in the framework of factorized scattering developed in [1] (for a later review see 
[2]). It was shown in [3] that in such theories using the scattering amplitudes as input it is 
possible to obtain a set of equations satisfied by the form factors. The complete system of 
form factor equations, which provides the basis for a programmatic approach (the so-called 
form factor bootstrap) was proposed in [4]. For a detailed and thorough exposition of the 
subject we refer to [5] on bulk form factors; later this approach was also extended to form 
factors of boundary operators [6, 7]. 

Although the connection with the Lagrangian formulation of quantum field theory is 
rather indirect in the bootstrap approach, it is thought that the general solution of the 
form factor axioms determines the complete local operator algebra of the theory. This 
expectation was confirmed in many cases by explicit comparison of the space of solutions 
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to the spectrum of local operators as described by the ultraviolet limiting conformal field 
theory [8, 9, 10, 11, 12]; the mathematical foundation is provided by the local commuta- 
tivity theorem stating that operators specified by solutions of the form factor bootstrap 
are mutually local [5]. Another important piece of information comes from correlation 
functions: using form factors, a spectral representation for the correlation functions can 
be built which provides a large distance expansion [13], while the Lagrangian or perturbed 
conformal field theory formulation allows one to obtain a short-distance expansion, which 
can then be compared provided there is an overlap between their regimes of validity [14]. 
Other evidence for the correspondence between the field theory and the solutions of the 
form factor bootstrap results from evaluating sum rules like Zamolodchikov's c-theorem 
[15, 16, 17, 18] or the A-theorem [19], both of which can be used to express conformal data 
as spectral sums in terms of form factors. Direct comparisons with multi-particle matrix 
elements are not so readily available, except for perturbative or 1/A^ calculations in some 
simple cases [3]. 

Therefore, part of the motivation is to provide non-perturbative evaluation of form 
factors from the Hamiltonian formulation, which then allows for a direct comparison with 
solutions of the form factor axioms. Another goal is to have a better understanding of 
finite size effects in the case of matrix elements of local operators, and to contribute to the 
investigation of finite volume [20, 21, 22] (and also finite temperature [23]) form factors 
and correlation functions. 

This program has been successfully pursued in the case of diagonal scattering theories 
(those without particle mass degeneracies), both in the bulk and with boundary [24, 25, 
26, 27]. However, an extension to the non-diagonal case is still missing, and the present 
work is a step in that direction. It is important to realize that non-diagonal theories, whose 
spectra contain some nontrivial particle multiplets (typically organized into representations 
of some group symmetry), such as the 0(3) nonlinear sigma model are very important for 
condensed matter applications (e.g. to spin chains; cf. [28]). The finite volume description 
of form factors can be used to develop a low-temperature and large-distance expansion for 
finite-temperature correlation functions [25, 29], which could in turn be used to explain 
experiments such as data from inelastic neutron scattering [30]. Therefore the extension 
to non-diagonal theories is an interesting direction. 

Sine-Gordon model can be considered as the prototype of a non-diagonal scattering 
theory, and it has the advantage that its finite volume spectra and form factors can be 
studied numerically using the truncated conformal space approach, originally developed 
by Yurov and Zamolodchikov for the scaling Lee- Yang model, but later extended to the 
sine-Gordon theory [31]. Its exact form factors are also known in full generality (at least 
in principle - see the discussion on multi-soliton form factors later), and so it is a useful 
playground to test our theoretical ideas on finite volume form factors. This cuts in another 
way too since our approach also provides a way to test the conjectured exact form factors 
in much more details than allowed by the usual methods listed above. 

The outline of the paper is as follows. We collect the necessary facts about the sine- 
Gordon form factors in section 2. In Section 3 we turn our attention to the breather form 
factors, which can be treated by the methods developed for diagonal scattering theories. 
Section 4 contains our results on soliton form factors, as well as a conjecture how the finite 
volume form factor formulae derived earlier in [24, 25] can be extended to non-diagonal 
theories. At present we are only able to perform a partial check of this conjecture due 
to difficulties in evaluating multi-soliton form factors numerically. Section 5 sums up our 
conclusions. The paper also contains an appendix devoted to a brief description of the 
algebraic Bethe Ansatz technique that can be used to obtain the finite volume description 
of multi-soliton energy levels. 
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2 Brief review of sine-Gordon form factors 



2.1 Action and S matrix 

The classical action of the theory is 

r /I rn^ 
A = J d^x i-d^^d^"^ + cos 

The fundamental excitations are a doublet of soliton/antisoliton of mass M. Their exact 
S matrix can be written as [1] 



where 



5^(^,0 = ^^(^,0^0(^,0 (2.1) 

sXz{e,i) = szX{e,o = ST{e,i) 



and 





dt sinh 2 
T sinh ^ cosh f 



J2sinh ^ilzilie— + (e""-?* - l) {e^i^-W2 + Q-n{i+0t/2^^ 



2 sinh 2|* cosh f 

(the latter representation is valid for any value of n G N and makes the integral represen- 
tation converge faster and further away from the real 6 axis). 

Besides the solitons, the spectrum of theory contains also breathers Br, with masses 

mr = 2Msm^ (2.2) 

The breather-soliton and breather-breather S'-matrices are also known, here we only quote 
the Bi — Bi amplitude that is needed in the sequel: 

sinhg + isin 

OBiBi(t^j = — — {2.6) 

smh u — t sm vr^ 

Another representation of the theory is as a free massless boson conformal field theory 
(CFT) perturbed by a relevant operator. The Hamiltonian can be written as 

H = J dx^: (dt^f + {d^^f ■.+fxjdx: cos/3<l> : (2.4) 

where the semicolon denotes normal ordering in terms of the modes of the /i = massless 
field. In this case, due to anomalous dimension of the normal ordered cosine operator, the 
coupling constant has dimension 

H ~ [massj ' 
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2.2 Breather form factors 



We consider only exponentials of the bosonic field $. To obtain matrix elements containing 
the first breather, one can analytically continue the form factors of sinh-Gordon theory 
obtained in [8] to imaginary values of the couplings. The result is 

FiV.i(^i,...,^n) = (0|e^"^*(°)|i?i(^i)...i?i(^„)> 



where 



Qa ^ ("^l 1 ■ ■ ■ 1 -^n) 

A(0 



i<j 

det [a + i- i]^ (T^^ljixi, ^^n) if n > 2 

sin7r,^a 



(2.5) 



sin 7T^ 



2 cos — Y 2 sin — exp 



dt t 



271 sin t 



and 



m 



v{i7r + e, -l)v{iTT + 9, -e)v{m + 6,l + e) 
vi-iiT - e, -l)v{-i7r - 6*, -^)v{-i7r -6,1 + ^) 

C 



^ ' 9 + i7r{2k + C)^ 



n 



L \e + iTT{2k - c) 



+ {N + 1- iVe~2*) 



exp 

2Nt+— siiih (^t 



iCe 



t \ 4 sinh I 27r cosh | 



2 sinh^ t 



(2.6) 



gives the minimal BiBi form factor^, while cr^"^ denotes the elementary symmetric poly- 
nomial of n variables and order k defined by 

n n 

n(a: + a:,) = ^x"-Vf (a;i,...,xj 

i=l 

Furthermore 



fc=0 



47r 
87r-^2 



or JVl. 



~1 477 




exp < 


[Jo T 



Air 



-2t 



+ 



sinh^ [fj] 



2 sinh 



t cosh 1 - 



t ) sinh t 



(2.7) 



is the exact vacuum expectation value of the exponential field [32], with M denoting the 
soliton mass related to the coupling /i defined in via [33] 

2-2A 



2r(A) / v^r (^) M 



vrr(l - A) I 2T (^) 



A 



57r 



(2. 



Formula (2.5) also coincides with the result given in [34]. Form factors of higher breathers 
can be constructed by considering them as bound states of Bi particles. Detailed formulae 
can be found in Appendix A of [35]. 



-•^The formula for the function v is in fact independent of A^; choosing N large extends the width of the 
strip where the integral converges and also speeds up convergence. 
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2.3 Soliton form factors 



At present, there are three independent constructions of sohtonic form factors available: 
the earliest one by Smirnov (reviewed in [36]), the free field representation by Lukyanov 
[37, 34] and the work by Babujian et al. [38, 39]. Here we use formulae from Lukyanov's 
work [34] . The simplest form factor for an exponential operator is with a soliton-antisoliton 
state and it takes the form 

(0 le-"*!""! S±(92)St(9i)> = - 

f|,(9) = ft(^)Gfflcot(f) "'-''(-P;^f (2.9) 

V^y ^sinhf^j 

where S+ and S- denote the soliton and antisoliton, respectively. The function G is given 
by the integral representation 



dt . .ff \ sinh(t(^-l)) 



-J exp <.j — sinh 




sinh (2t) cosh (t) sinh (t^) 
' si^h (0 n 9(0, ^^kf exp | J^^ ^e"^^* {l + N - N e"^*) 



xsinhM(l-!^V^ smh(t(e-l)) 



n J J sinh (2t) cosh (t) sinh (t^) 

p ^ (2fc+l+g)7r-i6) j p ^ 2k+i y Y ^ (2fc+l)7r-i6> j 



X- 



Y I {2k-l)n+ie \ Y f 2k~l+A ^ Y /' (2fc-l+g)7r+ie 



r r (2fc+2)7r-t6> \ -p / 2kTT+ie 



where, again, the second formula is eventually independent of the natural number A^; it 
provides a representation which converges faster numerically and is valid further away from 
the real 6 axis with increasing N. 

Higher form factors in all the available constructions are given in some quite complicated 
integral representation which we do not write down here. We have been unable to find a 
numerical evaluation of these analytic formulae at present, which precludes their use in a 
comparison with TCSA data. 



3 Breather form factors in finite volume 
3.1 A review of the theoretical predictions 

As breathers are singlet and so scatter diagonally, the formulae presented in the papers 
[24, 25] are directly applicable. The first ingredient is to describe the multi-breather energy 
levels corresponding to the states 

\Br,{0l)---BrAON)) 

whose finite volume counter part we are going to label 

ri...rjv,L 
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where the Ik are the momentum quantum numbers. In a finite volume L, momentum 
quantization is governed (up to corrections decreasing exponentially with L) by the Bethe- 
Yang equations: 

Qk{ei,..., On) = mr,L sinh Ou + Y^ ^r,r, {Ok " Oj) = 2nlk h e Z (3.1) 
where the phase-shift is defined as 

In practical calculations, because breathers of the same species satisfy an effective exclusion 
rule due to 

Sb^bM = -1 

it is best to redefine phase-shifts corresponding to them by extracting a — sign: 

so that Ssr{0) = can be taken for all s,r and all the phase-shifts can be defined as 
continuous functions over the whole real 6 axis. This entails shifting appropriate quantum 
numbers Ik to half-integer values. Given a solution 6i, . . . ,6fyr to the quantization relations 
(3.1) the energy and the momentum of the state can be written as 

AT 

E = nirf. cosh 9k 

k=l 

P = ^ rrir^ sinh 9k = — '^h 

k=l k 

(using that - with our choice of the phase-shift functions - unitarity entails Ssri9) + 
^rs{—9) = 0). The rapidity-space density of n-particle states can be calculated as 

Pr,...rAOi,---,0^)=detJ^-^ , jin)^ dQk{9u^..,9^) (3,2) 

The matrix elements of local operators between finite volume multi-particle states can be 
written as [24] 

I si...sm{{I'i, • • • , 0)l{A, • • • , lN})ri...rN,L\ = 



F?,,...s,r,...r^{hl + l7T,...,9[+tn,9,,..., 9^) 



\J Pri...rjv(^l, • • • , 9n)Psi...sm{9'i, . . . , 9'j^j) 



0(e-^'^) (3.3) 



(note that we cannot specify the phase of the matrix elements as it depends on phase 
conventions which can be different in finite and infinite volume) which is valid provided 
there are no disconnected terms. Such terms appear when there are two breathers of the 
same species in the two states whose rapidities coincide. Apart from some very special 
cases, this happens only when the two states are identical, which is called the diagonal 
matrix element. The necessary formulae are written down in [25]; instead of quoting here 
the rather lengthy general expression, we present the two particular cases we need: 

i({J}|O(0,0)|{/})i,i = J^%^ + ((9) + 0(e-'^'^) (3.4) 
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({/l,/2}|O(0,0)|{Ji,J2})ll,L 



miL cosh 6' 

^n(^i,^2) (cosh^i + cosh^2) i^n(«7r,0) 



Pll(^l,^2) 

+ {O) + 0(e-^'^) 



6 



where 

(O) 

is the vacuum expectation value of the local operator O and 

J'liiOi, 62) = limFi^ii(^2 + 111 + 6,61+111 + e, ^1, 62) 

is the so-called symmetric evaluation of the four-particle form factor. Note that in the di- 
agonal case there is no possible phase difference between the finite- volume matrix elements 
and the infinite volume form factors, as any phase redefinition of the state drops out from 
the matrix element. 

The above predictions for finite volume energy levels and matrix elements are expected 
to be exact to all (finite) orders in 1/L [24, 25] (note that the exponential corrections are 
non-analytic in this variable). 

3.2 Numerical results 

To evaluate the form factors numerically, we use the truncated conformal space approach 
(TCSA) pioneered by Yurov and Zamolodchikov [40]. The extension to the sine-Gordon 
model was developed in [31] and has found numerous applications since then. The Hilbert 
space can be split by the eigenvalues of the topological charge Q (or winding number) and 
the spatial momentum P, where the eigenvalues of the latter are of the form 

27rs 

s is called the 'conformal spin'. 

In sectors with vanishing topological charge, we can make use of the symmetry of the 
Hamiltonian under 

C: $(x,t) ^ -$(x,t) 

which is equivalent to conjugation of the solitonic charge. The truncated space can be split 
into C-even and C-odd subspaces that have roughly equal dimensions, which speeds up the 
diagonalization of the Hamiltonian by roughly a factor of eight (the required machine time 
scales approximately with the third power of matrix size). We remark that there is another 
discrete symmetry 

V : $(a;,t) -$(-x,t) 

corresponding to spatial parity, which is a symmetry in s = sectors. However, we use the 
s 7^ states extensively in our calculations. 
For the calculations, we choose the operator 

O =: e*^* : 

which is essentially one half of the interaction term in the Hamiltonian in (2.4). The 
semicolons denote normal ordering with respect to the A = free massless boson modes. 
This operator has conformal dimension 

Svr 

Using relation (2.8) we can express all energy levels and matrix elements in units of (ap- 
propriate powers of) the soliton mass M, and we also introduce the dimensionless volume 
variable / = ML. The general procedure is the same as in [24, 25]: the particle content 
of energy levels can be identified by matching the numerical TCSA spectrum against the 
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Figure 3.1: One-particle form factors at ,^ = 50/239 



predictions of the Bethe-Yang equations (3.1). After identification, one can compare the 
appropriate matrix elements to the theoretical values given by (3.3) and (3.4). 

Due to level crossings, at certain values of the volume L there can be more than one 
TCSA candidate levels for a given Bethe-Yang solution; identification can be completed 
by selecting the candidate on the basis of one of the form factor measurements, which 
still leaves other matrix elements involving the state as cross-checks. Level crossings also 
present a problem in numerical stability, since in their vicinity the state of interest is 
nearly degenerate to another one. Since the truncation effect can be considered as an 
additional perturbing operator, the level crossings are eventually lifted. However, such a 
near- degeneracy greatly magnifies truncation effects on the eigenvectors and therefore the 
matrix elements [27]. This is the reason behind the fact that in some of the figures there 
are some individual numerical points that are clearly scattered away from their expected 
place. 

The simplest matrix element to test is the theoretical prediction (2.7) for the exact 
vacuum expectation value, and we did check it against the numerical results. However, it 
has already been subjected to extensive verification against TCSA in [41], therefore we do 
not dwell on this issue here. 

One-particle form factors can be tested using a reorganized form of the relation (3.3): 



(0|e*^*W|fi„(0))| = v/;^(0|e*^*W|{0}) 



Ofe-'^'^) 



(the form factors themselves are independent of the rapidity of the particle due to Lorentz 
invariance). For the lowest three breathers the agreement between numerics and theory is 
illustrated in figure 3.1. 

The next interesting matrix element is the two-particle one, for which we can use the 
relation 



F.'iie^J,) = (0|e*^*(°)|i?i(e~i)i?i( 



-0(e 



which is tested in figure 3.2. This provides a test of the two-particle form factor for real 
rapidity differences. However, the two-particle form factor function also appears in the 
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(off- diagonal and diagonal) Bi-Bi matrix elements 



,({/}|e^^*(°)|{/}>, 



m\L\/ cosh Q' cosh Q 



0(e 



F<?(zvr,0) 



miLcosh6' 



~+e^i(/3) + 0(e"^'^) 



which makes it possible to test the form factor for rapidity differences with imaginary part 
vr as shown in figures 3.3 and 3.4. We remark that in plots against rapidity the small- 
rapidity deviations are due to truncation errors, while the ones for large rapidity result 
from the neglected exponential corrections. 

We can also test the three and four particle form-factors using By — ByBx and off- 
diagonal -Bi-Bi — -Bi-Bi matrix elements, as shown in figures 3.5 and 3.6. Diagonal -Bi-Bi — 
BxBx matrix elements reveal some complications, and are treated in subsection 3.4. 

All of these tests show excellent agreement between numerics and theory, which make 
us confident both in the finite volume form factor formalism and the correctness of the 
exact form factors (2.5) predicted by the bootstrap. 



3.3 B2 matrix elements and ^u-terms 

We can also use the identified B^ state to compute matrix elements involving it. However, 
in some cases simply using (3.3) gives a rather poor match. It turns out that accounting 
for some of the hitherto neglected exponential corrections can improve this; the leading 
corrections come from so-called /i-terms which were constructed in [42] building upon the 
ideas in Liischer's seminal paper [43]. We remark that the method outlined below eventually 
does more; it entails a resummation of these correction beyond the mere construction of 
the leading exponential correction. 

Using the ideas in [42], we can model a finite volume B2 state as a pair of B\ particles 
with complex conjugate rapidities by solving the B\B\ Bethe-Yang equations (here written 
in exponential form): 

gimiLsinh(e±m)^^^^^(^2m) = 1 (3.5) 

The solution for u has the large volume behaviour 

«^ !l + Ce-'^^™^^^ (3.6) 

with 

\~ mf . vr^ 

/i = y m\ = mi sm — 

V 4 2 

and C some numerical constant. Therefore, for L — )■ oo one obtains the usual bootstrap 

identification 

B^iQ) ^B,{e + m^/2)B,{e - 
and the momentum and energy of the state tends to 

mi sinh(^ + i-u) + mi sinh(^ — 2m) = 2mi cosM sinh^ — )■ m2sinh^ 

L— >oo 

mi cosh(6' + in) -|- mi cosh(6' — m) = 2mi cosncosh^ — )■ m2Cosh6' 

L— >-oo 

where 

m2 = 2mi cos — 
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Figure 3.2: BiBi two-particle form factors. The first plot shows spin-0 data with ^ 
50/239, while the second one shows spin-1 data with ^ = 2/7. 
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Figure 3.5: Bi-BiBi matrix elements at ,^ = 50/239. 
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Figure 3.6: Off-diagonal BiBi-BiBi matrix elements at ^ = 50/239. The individual 
deviations are caused by the proximity of level crossings. 
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is the mass of B2 in terms of that of Bi, consistent with (2.2). 

In order to do this, however, formula (3.3) must be continued to complex rapidities, 
which requires matching the phases of the matrix elements on the two sides. This is not 
difficult to do by observing that the TCSA matrix elements are all reaP, while the phase 
of the minimal form factor function f^{6) (2.6) is just half the phase of the two-particle 
5- matrix SbiBi{(^)- Using the results of [42] one obtains 



(010(0, 0)|{/}),. = ^ V^Ss.sA -2^u)FOie + ^u,e-^u) ^ 

Y Pii{^ + iu,9 - iu) 

where 9, u is the solution of (3.5) with the correct momentum, ie. 

miL sinh(6' + iu) + rriiL sinh(^ — iu) = 2iil 

and the dots denote further (and generally much smaller) exponential corrections; the ± 
sign accounts for the remaining phase ambiguity from the square root. One can similarly 
evaluate Bi — B2 matrix elements using 



l({/'}|0(0,0)|{/}) 



v/5B,B,(-2zn)F5i(z7r + 6', 6 + zn, 6 



lU] 



2,L 



rriiL cosh^' pii{6 + iu,9 — iu) 



where the absolute values are necessary because we have not compensated for the phase 
difference due to the presence of the additional B^. We observed substantial improvement 
over the naive matrix elements (3.3) that neglect the yU-term contributions; some of our 
data are demonstrated in figure 3.7. The effect is very dramatic for the Bi — B2 case, the 
naive prediction drastically disagrees with the numerical results, while the inclusion of the 
/i-terms leads to excellent agreement. This is in fact a bit unexpected, as the value for the 
exponent for ^ = 2/7 is 



777 

fiL = dmj- 0.376511 , I = ML 

whose coefficient is not particularly small. In the original example in [42] one of the 
particles was loosely bound, resulting in a very small value for fi, which in turn led to 
large finite size corrections from the particle splitting up in finite volume. This is not case 
for B2 considered as a bound state of two BiS at the particular value of the coupling in 
question. However, the full behaviour is determined not only by the exponent, but also by 
the detailed behaviour of the ^-matrix amplitude and the form factor near the pole; it is 
the latter that enhances the correction in the Bi — B2 case. We return to this observation 
and its implications in the conclusion. 



3.4 Diagonal BiBi — BiBi matrix elements 

For the diagonal matrix element we observe significant difference between the numerical 
results and the theoretical predictions as shown in figure 3.8. It is apparent from the data 
that agreement is better for lower values of ^. For large volumes, this is easily explained 
by the better convergence of TCSA since the dominant source of errors in the large volume 

^This is the case because the TCSA representation of the Hamiltonian turns out to be a real and 
symmetric matrix, due to the fact that the matrix elements of operators e''^'^* are all real. As a results, 
all eigenvectors are also real and therefore all the matrix elements computed from TCSA are real as well. 
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Figure 3.7: Naive predictions for matrix elements involving B2, together with ones including 
/i-term corrections compared to numerical data for the i?2-vacuum and B2 — Bi matrix 
elements at ^ = 2/7. 

data is the Hilbert space truncation, and improvement of truncation errors with decreasing 
^ was observed in all the data we analyzed. 

However, this cannot be the origin the deviation in the medium range of volume (5 < 
/ < 15). This deviation also decreases for smaller ^, and is likely to have its origin in a 
/i-term correction. Evaluating such a /i-term requires modeling Bi as a soliton-antisoliton 
bound state, which we are not able to do at present as it necessitates numerically handling 
eight-particle solitonic form factors. As for the variation of the /t-term with ^, analogous 
behaviour was observed for the B2 case, where we also have analytic confirmation of the 
numerical results. 

3.5 Higher breathers 

We also tested numerous other form factors including higher breathers (among them B2 — 
BiBi, B2-B2, Bi — B3, B2 — -B3 and B3 — B3 matrix elements), and in all cases found good 
agreement between theoretical predictions and TCSA data. These tests cover all multi--Bi 
form factors (2.5) up to 6 particles (taking into account that Bn is obtained as a fusion of 
n first breathers, as described in Appendix A of [35]). 

4 Soliton form factors in finite volume 

4.1 Finite volume form factors in non-diagonal theories 

We can follow the lines of the reasoning in [24] that lead to the formula for the finite volume 
form factors. We recall the main outlines of the arguments, making modifications where 
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Figure 3.8: Diagonal -Bi-Bi — -Bi-Bi matrix elements for ^ = 2/7 and ^ = 50/311. 
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necessary. We consider the spectral representation of the Euchdean two-point function 

(o(x)o'(o,o)) = EE n/ If c„(^i.^2,...,^„)x 

i^f„,jei,^2,...,^n)+exp (^Rjyi^^coshe^ (4.1) 

where 

i^f .^J^i, ^2, . . . , ^n)+ = n....„(^i, • • • , 0^\O\Q, 0)|0) = i^f.'..,„(^i + zvr, ^2 + ZTT, . . . , ^„ + ztt) 



(which is just the complex conjugate of F^' for unitary theories) and R = \/t^ + x"^ is the 
length of the Euclidean separation vector x = (r, x). 

In finite volume L, the space of states can still be labeled by multi-particle states but 
the momenta (and therefore the rapidities) are quantized. Denoting the quantum numbers 
/i, . . . , J„ the two-point function of the same local operator can be written as 

oo 

{O{t,0)O'{0,0)), = EE E(0|^(0'0)l{^i'^2,...,/n})?x (4.2) 

n=0 r h...In 



(r 



•) {{h, h,..., In}\O'{0, 0)\0)l exp -r E ^i, cosh 6, 



k=l 



where the index r enumerates the eigenvector of the n-particle transfer matrix (as an exam- 
ple cf. Appendix A where these are constructed for sine-Gordon solitons); these are usually 
rapidity dependent linear combinations in the space of particle multiplet indices ii . . .in- 
We can assume that the wave function amplitudes (A.l) of these states are normalized and 
form a complete basis: 

E<U({^a)^i:U(R})*=5.. 

r 

In writing (4.2) we also supposed that the finite volume multi-particle states I2, ■ ■ ■ , In})r,L 
are orthonormal and for simplicity restricted the formula to separation in Euclidean time 
r only. The index L signals that the matrix element is evaluated in finite volume L. Using 
the finite volume expansion developed by Liischer in [43] one can easily see that 

(0(r, O)O'(O)) - (0(r, 0)O'{0))l ~ ©(e"^^) (4.3) 

where n is some characteristic mass scale. 

We can now rewrite the infinite volume correlator in the basis of multi-particle transfer 
matrix eigenstates 

°° / " roo jQ \ 

{O{r,0)O\0M = EE n/ ^]FA0i,...,0nrx 

F^' {9^,..., 9nY^^-^ exp (-r E cosh 9^) 

where 
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The remainder of the argument follows the lines of the paper [24]. Essentially, we compare 
the discrete sum with the integral, and realize that up to exponentially small terms in L 
the relation is given by the Jacobian of the mapping between the quantum number and 
the rapidity space, i.e. the density of states. We obtain 



(o|o(o,o)i{/i,...,/j)i:') 



Ofe^^'^l 



(4.4) 



where^ 



is the density of states with internal wave vector as given in (A. 4), and 9^ are the 
solutions of the Bethe-Yang equations (A. 3) corresponding to the state with the specified 
quantum numbers Ji, . . . , J„ at the given volume L. 

This can be easily generalized to matrix elements with no disconnected pieces: 



(^) 



({/(,..., /;,}|0(o,o)|{/i,...,M) 



(r) 



?0( 



hi, ■ ■ ■ 1 0i\9i, . . . , On)^"''^^ 



+ 0(e-'^'^) 



(4.5) 



where 



F^(^~;„...,^~;i^~i,...,^~^r)(^'^) 



where the bar denotes the antiparticle. 

This can be easily applied to soliton-antisoliton states. The algebraic Bethe Ansatz 
gives two eigenvectors in this subspace (cf. A. 3): 



5(Ai)fi oc 




for Ai 
for Ai 



h6»2 



h6»2 



*2 



I.e. 



vI/W = -^(0,+l,+l,0) 

^(-^ = ^(0, +1,-1,0) 

Since the value of the magnonic variable Ai is explicitly known, it can be eliminated from 
the Bethe-Yang equations, resulting in the following quantization conditions for the soliton- 
antisoliton pair: 



Q[^\dud2) = MLsmhOi - zlog5±(^i - 62) = 27r/i 
<5?^(^i,^2) = MLsinh^2 -«log>S±(^2 -^1) = SvrJa 



(4.6) 



^We note that 11' is not necessarily the same scale as in (4.3) and its value depends on the spectrum 
and bound state fusion angles of the model. 
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where S± are the eigenvalues of the two-particle 5*- matrix given in (A. 2) and the ± dis- 
tinguishes the two states. It is also easy to eliminate the magnon Ai from the density of 
states to obtain 



(901 
801 



002 
802 



From (4.4) we obtain 



(0|0(0,0)|{/„/2})l^) 



0(e 



-HL\ 



(4.7) 



where 



in terms of (2.9) and 9i^2 are the solutions of (4.6) at the given volume L with quantum 
numbers Ji 2- 

At the moment, there is no general theoretical result for matrix elements with discon- 
nected pieces (e.g. diagonal ones). However, for one-soliton matrix elements one can easily 
write down the appropriate generalization of eqn. (3.4): 



.({/}|0(O,O)|{J})_ 



FA-m) 



ML cosh 9 

where F_+ is given in (2.9) and the rapidity 9 is quantized as 

ML sinh^ = 2txI 



(4. 



4.2 Numerical verification of solitonic matrix elements 

The formulae derived in the previous subsection show an excellent agreement with TCSA. 
Figure 4.1 demonstrates this for the matrix element between the vacuum and the (sym- 
metric or antisymmetric) soliton-antisoliton two-particle states as given in eqn. (4.7). We 
can also study one-soliton-one-soliton matrix elements. For the non-diagonal ones we can 
use the = M = 1 case of (4.5) while for the diagonal ones it turns out we must slightly 
modify (4.8) by changing the sign of the F_+ term: 

+({/}|o(o, o)|{/})+ = - f;;^"7l + gm (4.9) 

ML cosh 9 

It is not yet clear to us what is the reason behind this, but, as shown in figure 4.2 this 
produces an excellent agreement with the TCSA data. It is possible that the bootstrap 
solution (2.9) is off by a sign^, or that the sign here is related in some way to the crossing 
of a soliton. 



5 Conclusions and outlook 

In this work we investigated the comparison between form factors of sine-Gordon theory 
obtained from the bootstrap and finite volume matrix elements given by the truncated 

^This would not be detected by checking off-diagonal matrix elements, since due to phase differences 
these compare only the absolute values according to (4.5). 
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Figure 4.1: Vacuum-sohton-antisohton form factors at ^ = 2/7. The individual deviations 
are caused by the proximity of level crossings. 
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Figure 4.2: Soliton-soliton form factors: off-diagonal and diagonal case at ^ = 50/239. 
The cusp in the continuous line is an artifact caused by taking the absolute value. 

conformal space approach. Generally speaking, we found an excellent agreement between 
the theoretical predictions and the numerical results. We also generalized the formalism 
developed in [24, 25] to non-diagonal scattering theory in order to include solitonic matrix 
elements in the test. 

In this latter respect, however, several problems are still open. One of them is to 
find a suitable formula for matrix elements with disconnected contributions in the case 
of non-diagonal scattering (at present we only have the simplest, one-particle-one-particle 
case covered). The other, more technical problem is to find a suitable framework for the 
numerical evaluation of the multi-soliton form factors; until then, the full strength of the 
algebraic Bethe Ansatz formalism presented in Appendix A cannot be utilized. In addition, 
the solution of the problem in the case of diagonal BiBi — BiBi matrix elements, described 
in subsection 3.4, also hinges on the ability of handling such matrix elements. A further 
issue is the presence of the minus sign in (4.9) which requires a detailed look into the soliton 
form factor bootstrap, which is out of the scope of present work. 

Another interesting issue is the role played by the so-called /i-terms, which are expo- 
nential corrections to the volume dependence of the matrix elements resulting from particle 
fusions (in other words, three-particle couplings between on-shell particle states). The ex- 
ponential decay of these corrections depends only on the mass spectrum; however, their 
strength is also determined by appropriate form factors, as illustrated in subsection 3.3 by 
the fact that for the same B2 state (hence the same /i-term kinematics), the significance 
of the corrections depends very much on the matrix element the state is involved in. 

In this work we generally took the stance of using the bootstrap results to predict the 
numerical finite volume results. However, one could take the opposite route and extract 
field theoretically interesting matrix elements from the finite volume data (which we also 
did for cases involving two-particle matrix elements). This is common in lattice field theory 
and was advocated by Lellouch and Liischer in [44] as a way to circumvent the so-called 
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Maiani- Testa no-go theorem [45]. We have previously used a similar approach to determine 
resonance parameters from finite volume spectra [46]. 

The /i-term results have an implication for such studies, which aim to determine matrix 
elements (e.g. related to weak interaction) from finite volume data (whether using lattice 
Monte-Carlo or some other numerical method). The important lesson is that even if the 
appropriate exponential factor fi (as determined by the mass spectrum) is not small and 
therefore volume suppression is expected to be strong, the /x-term can still have a sizable 
effect at the values of volume the applied numerical method (whether TCSA or lattice) 
can reach. In addition, this term could vary strongly between different matrix elements 
involving the same state, thereby producing very large systematic errors when extracting 
matrix elements of local operators from numerical simulations (indeed, the corrections 
could of the order of 100% as illustrated by the B2 — Bi matrix element plotted in figure 
3.7). 
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A Algebraic Bethe Ansatz for multi-soliton states 

A.l The sine- Gordon soliton Bethe- Yang equations and the multi- 
particle transfer matrix 

Our starting object is the A^-particle monodromy matrix that can be written as 
or as a 2 X 2 matrix in the a, h indices 



M (A| {^1, . . ..e^}t:!l;X = '^Sf (A - 0i)szi:ix - 02) ■ • •5^,.,(A - e, 



M{x\{e,}) 



A{X\{9,}) B{\\{9,}) 
CiX\{9,}) Di\\{9,}) 

The operators A, B, C and D act in the 2^-dimensional "isospin" space (the space spanned 
by the soliton doublets) 



n 

2 



V;v = (g)C 



k=l 



which is spanned by the basis 
One can introduce the vector 



corresponding to all solitons being positively charged. The general transfer matrix 

r (A| {9,}) = M (A| {9,}): = A (A| {9,}) + D (A| {9,}) 

gives rise to the specified transfer matrices 

rAm) = r{X = 9,\{9,}) 

In terms of these the ra-particle quantization relations, ie. the Bethe- Yang equations on 
the circle can be written as [47] 

where \E' {{9k}) € Vat is the wave-function amplitude vector. 
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A. 2 The Algebraic Bethe Ansatz 

We use the well-established machinery of the algebraic Bethe Ansatz (for a pedagogical 
introduction see [48]). The sine-Gordon S'-matrix satisfies the Yang-Baxter equations 

where 6ij = 9i — 6j. As a result, the monodromy matrix satisfies the relations 

- (A| {e,}t\ M (/.I {e,}ti = M (/ii {e,})i m (a| {6,})% siit^{\ - /x) 

Therefore the transfer matrices form a commuting family 

T (A| {9,}) T (/i| R}) = T (/i| R}) T (A| {9,}) 

and can be diagonalized simultaneously with common eigenvectors for all A: 

r (A| {9u}) ^iP^ {{9u}) = A JA| {^4) {{9k}) 

where Aq (A| {9k}) G are the eigenvalues, ipa {{9k}) G Vn are the eigenvectors and 
a = 1 ... 2^. The independent solutions of the Bethe equations for the wave vector are 
then given by 

* {{Ok}) = i^a {{Ok}) 

provided the particle rapidities satisfy the Bethe- Yang equations 



e 



iMLsinhl 



'^A,(^,|{^fc}) = l J = 1...N 




In this way we reduced the problem to finding the eigenvalues of the transfer matrix and 
then solving a system of coupled scalar equations. 

Note that the sine-Gordon scattering preserves the solitonic charge Q and charge parity 
C, which act on Vn as follows: 

Qvi^...iN = 

Cvi,,„i^ = Vi^,„i^ where i = -i 

The transfer matrix T commutes with these operators. In addition, B and C act as 
lowering and raising operators in charge space 

[Q,B{X\{9k})] = -2i?(A|{^,}) 
[Q,C{\\{9k})] = +2C{\\{9k}) 

We now look for the eigenvectors of the transfer matrix in the form 

* ({AJ I {Ok}) =U^B (All {9k}) ...B{K\ {9k}) n (A.l) 

where A^ are the so-called 'magnons' and A/^i^ is some normalization factor. They satisfy 
the eigenvalue equations 

r{\\{9^,..., 9„}) VP {{K} I {9k}) = A (A, {AJ | {9^}) * ({AJ | {9^}) 

provided the algebraic Bethe Ansatz (ABA) equations hold 



^ n(\ \ \ 1 sinh ^ 

n «(^. - = n ^^^1 = = 

t=i fc^^«(Afc-Aj) ^r(t^,4) sinh (I 
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The eigenvalue is given by 

/ r N 

A (A, {A J I {^1,...,^^}) = [l[ST{\k- \)-' + l[ST{\-ek,ollsT{\- \k)-' 

.k=l 
N 

l[Soix-ek,o 



Kk=l k=l k=l 

N 

X 

k=l 



Notice that a{6) (and therefore the ABA equations) has the natural periodicity 

a{e + zttO = a{e) 

A. 3 The two-particle case (A^^ = 2) 

The only interesting eigenvectors are the ones in the Q = subspace, for which we can 
write the Ansatz 

5(Ai)fi 

The single ABA equation takes the form 



sinh fi^L^^±ii\ sinh f 
a(Ai-^i)a(Ai-^2) - V ^ / V ^ 



sinh ( sinh 



and (up to periodicity) has two independent solutions 

(1) : Ai = — - — + t- 

(2) : Ai = — - — + ^7r^— 
We can then evaluate the corresponding eigenvalues: 

A (A. {A,} I {e,M) = (^^(^ + '^"'^'srT-t)''''^ ) - - "'^ 

Putting A = (to get the Bethe-Yang phase-shift for the first particle) we find 

(1) : Ai9^,{\,}\{9^,92}) = -S+i9,-92) 

(2) : Ai9^,{X,}\{9^,92}) = -S49,-92) 

where 

S±{9) = {St{9)±Sr{9))So{9) (A.2) 
are the eigenvalues of the two-particle S matrix in the neutral subspace. 

A. 4 The four-particle case {N = 4) 

Albeit we are not able to use them for detailed numerical comparison yet, we briefly review 
the four-particle results to give a better understanding of what follows. Suppose that the 
four rapidities are ordered as 

9i<92<9s< 9i 
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A. 4.1 Q = 2 sector 

This subspace is four dimensional, and with the Ansatz 

B{Xi)n 

there are two types of solutions: 

• Al = yU + z| 

There are three such solutions, typically one with /x the between 62 and 6*3, one around 
61 and another one around 64. 

• Ai = /i + ^(l + Of 

There is a single solution, with n the between 62 and ^3. 



A.4.2 Q = sector 

There are 6 eigenvectors here, which can be classified by their parity under charge conju- 
gation C. The Ansatz is 



B{\i)B{\2)n 

and the ABA equations are 

fe=i fe=i 

For the three C = — 1 eigenvectors, the magnons take the form 

Ai = /ii + y A2 = /i2 + 

with the positions 

61 < jii < 62 < O'i < 112 < O4 
or 9i < ij,2 < O2 < 63 < jj,! < 64 
or 6*1 < 6*2 < /ii,Ai2 < 6*3 < 6*4 
There are three C = +1 eigenvectors. One has magnons of the form 



ITT ITT 

Ai = /ii + y A2 = /i2 + y 



and another one with 



with the typical positions 

9i < fii < 02 < 9-3 < IJ^2 < Oa 
The last eigenvector is of the type 

in{l + x) ivrfl — x) 

Ai = /i H Ai=/i + 



2 ' ^ 2 

where /i is exactly 

^1 + ^2 + ^'3 + 6'4 
^= 4 

and the fundamental range of x is 

< a; < ^/2 

given that the algebraic Bethe Ansatz has periodicity itt^ (and the sign of x does not 
matter). 
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A. 5 The density of states 

The quantization conditions for the multi-sohton states can be written as follows: 



e 



'^A(^„{Ai,...,A,}|{^i,...,^^}) = 1 , j = l,...,7V 



N 



k=l k=/=j 



a{\j - Afe) 



where is the number of solitonic particles and Q = — 2r is their total topologi- 
cal charge. The magnonic configuration as detailed in the previous subsection fixes the 
"isospin" structure of the state in the 2^ dimensional internal charge space. The two equa- 
tions must be solved simultaneously. In principle, the positions of the magnons are fixed 
in terms of the rapidities 6i, . . . ,6^ so the density of states with a given A^, r and fixed 
isospin structure can be computed from the quantization relations of the solitons 

QjiOi, . . . ^7v|Ai, ...,K) = ML sinh^j - i log A (^j, {Ai, . . . , AJ | {^i, . . . , On]) = 27r/j 

by taking the Jacobi determinant of the rapidity i— t- quantum numbers mapping [24] 

'dQj 



p(6'i, . . . ,6'Ar) = det 



j,k=l,...,N 



where it is understood that we differentiate also the dependence of the A^ with respect to the 
rapidities. However, there is an easier way to obtain the density of states by extending the 
set of BY equations by those of the magnons, considering also the magnonic quantization 
relations 



N r 

Mj{9i, . . . 9n\>^i, ■ ■ ■ , K) = log a{Xj - 9k) - log 

k=l k=/=j 



a(Afc - Xj) 
a{Xj - Xk) 



Then one has the following result 



di9,X) /d^^^ 



(A.4) 



Xs — ^As 



where the first determinant is the (N + r) x (N + r) Jacobi determinant formed by including 
the magnonic equations, differentiating by keeping the 9k and Xj independent, and the 
second is the r x r one containing only the magnonic relations. This can be proven with a 
bit of labour by expressing the derivatives 

aA, 

d9k 

by differentiating the magnonic relations and using simple properties of determinants. 
However, the result is very easy to understand intuitively. Since the magnons do not have 
independently chosen quantum numbers (their position is fixed by the rapidities of the 
solitons), the solitonic phase space volume can be obtained from the phase space volume 
calculated with all variables (rapidities and magnons) involved, divided by the volume 
contribution of the magnons themselves. 
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